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ABSTRACT 


C3 


~ The  four-pole  method  common  to  electrical  circuit  analysis  is 
applied  to  analysis  of  sound  transmission  through  layered  media.  This 
method  provides  a characterization  of  the  system  by  means  of  a 2 x 2 
matrix,  called  the  four-pole  matrix  of  the  system.  A general  four-pole 
matrix  is  derived  for  a single  viscoelastic  layer.  A sound  trans- 
mission system  consisting  of  many  layers  is  analyzed  by  specializing 
this  matrix  for  each  layer,  and  multiplying  the  resulting  matrices, 
giving  an  equivalent  four-pole  matrix  for  the  system.  Electrical 
analogs  are  easily  computed  from  the  terms  of  the  four-pole  matrices. 


BRIEF  DISCUSSION  OF  THE  FOUR-POLE  PARAMETER  METHOD 


\ 


The  problem  considered  is  the  determination  for  a linear  mechanical 


I , system  represented  in  Figure  1 of  the  input  force  and  velocity  (F-p  V-, ) 

in  terms  of  the  output  force  and  velocity  (F2,  V2)>  and  quantities 
that  depend  on  the  character  of  the  system.  If  steady-state  harmonic 
inputs  and  outputs  only  are  considered,  the  relationships 

F1  = a jjl  F2  + « 12  V2 
V1  = a 2i  F2  + a 22  V2 

are  valid'^,  where  the  or  i j depend  on  the  system  under  consideration. 
The  a ij  are  called  the  four-pole  parameters  of  the  system.  Determin- 
ation of  the  a ij  is  a relatively  simple  matter.  It  is  merely  assumed 
that  the  input  and  output  have  representations*  A e*-wt  where  A.,  <*> 
are  constants  (A  is  the  complex  amplitude  and  u>  the  angular  frequency), 
t is  time,  i is  vCT.  The  equations  for  velocity  and  force  are  written 
by  Newton's  laws  and  constitutive  relations  of  the  system,  and  evaluated 
at  input  and  output  terminals  of  the  system,  and  the  a ij  are  solved 
for  algebraically.  A simple  example  of  a spring  - mass  system  is 
presented  to  illustrate. 


For  the  system  represented  in  figure  2.,  Newton’s  second  law  yields: 


F1  - F2  - M 3? 


The  constitutive  relation  for  the  spring  is: 


F2  = K(X1  - X2) 


(2) 


* For  an  explanation  of  this  representation,  see  reference  (2)  sect.  3-3. 
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It  is  now  assumed  that  the  input  and  output  forces  and  velocities 

are  sinusoidal  with  constant  amplitude,  i.e.,  F-i  = A-t  e-^0*-  y.  = B-iCiwt 
j = 1,2.  Then  from  (1)  and  (2):  J J 3 J * 

F-^  - F2  = iumV-^ 
iwF2  = K(V1  - V2) 

Solving  for  Fj_,  V-j^  : 


F1  = (1 £r—  ) F2  + i°in  V2 

V1  = ~ f2  + v2 


Using  relations  (3)  and  (4),  or  their  inverse: 


\vl  = (/,iJ  } | v1!  ' (/JiJ  ) = <aiJ) 


-1 


(3) 

(4) 


any  two  of  the  four  quantities  may  be  computed  if  the  other  two  axe 
known.  For  example,  if  the  system  is  driven  as  shown  in  figure  3,  then 
V2  - 0,  Fj_  = A e:Uut . Then  from  (3)  and  (4): 

^ f f2  v2 


L 


F^  = A e 


F1  Vl 

iu>t 


Figure  3 


: 

( 


Referring  to  Figure  4,  the  force  voltage  analogy*  can  be  obtained 
from  the  four-pole  parameters  by: 


Figure  4 


the  use  of  the  impedance  formulae: 
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Mechanical 

Electrical 

Input  force  Fj_ 

Input  voltage  F-^ 

Output  force  F2 

Output  voltage  F2 

Input  velocity 

Current 

Output  velocity  V2 

Current  V2 

Table  1 


For  the  above  problem: 


iwm 


If  it  is  specified  that  V2  = 0,  FL  = A eiwt,  the  force-voltage 

analogous  circuit  for  the  system  shown  in  Figure  3 is  as  shown  in 
Figure  5.  nu 


Figjire  5 

For  n systems  connected  in  tandem,  the  kth  system  from  the 
input  terminal  being  characterized  by  the  four-pole  matrix  the 

equivalent  single  four-pole  matric  (arij)  for  the  entire  system  is*1): 

m 


arij 


k=l 


(8) 


5 


vl 


It  is  noted  that  the  order  of  multiplication  of  the  four-pole 
matrices  in  formula  (8)  is  required  to  be  the  same  as  the  order  of 
connection  of  the  subsystems,  going  from  input  to  output. 


DERIVATION  OF  THE  FOUR- POLE  PARAMETERS  FOR  A 
GENERAL  VISCOELASTIC  LAYER 

Consider  a layer  of  viscoelastic  material  acted  upon  by  sound 
pressure  fields  as  shown  in  Figure  7,  where 


T ' 


1 


Pi*  V1 


T 


11= 


p2,  v2 


Figure  7 


p represents  pressure,  and  v represents  velocity.  Using  a four-pole 
technique  similar  to  that  explained  in  the  previous  section,  we  can 
compute  a four- pole  matrix  ( aij)  such  that: 


= («ij) 


P2 
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i 
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Furthermore,  an  electrical  analog  for  the  layer  can  be  computed  using 
formulae  (5)*  (6),  and  (7).  Electrical  analogs  are  commonly  used  in 
acoustical  problems,  and  can  be  obtained  by  other  means;  see  for 
example  reference  (4).  The  most  important  advantage  of  using  a four- 
pole  technique  accrues  when  a plute  is  composed  of  layers,  which  may 
have  different  viscoelastic  properties.  The  equivalent  four-pole 
(and  subsequently  the  electrical  analog)  for  the  entire  composite  is 
obtained  rather  simply  by  application  of  formula  (8).  The  four-pole 
parameters  are  computed  below  for  a "typical"  layer,  and  the  four-pole 
matrix  for  a specific  layer  is  obtained  by  substituting  the  appropriate 
parameters  into  this  "typical"  matrix.  It  is  emphasized  that  the 
analysis  that  follows  is  only  valid  for  plane  waves  of  pure  harmonic 
character  propagating  through  linear  viscoelastic  materials  normal 
to  the  plane  of  the  layer. 

The  equation  of  motion  (9)  and  the  constitutive  relation  (10) 
under  the  assumptions  above  for  the  layer  in  Figure  7 are: 


X = the  lagrangian  coordinate  normal  to  the  plane  of  the  plate 
u = displacement  in  the  X direction 
M = a complex  modulus*,  M =1^(1  + ir|) 

M = the  appropriate  elastic  modulus  of  the  system 

i = v^r 


(9) 

(10) 


h = a loss  factor  depending  on  the  material,  its  temperature,  and 
the  frequency  of  excitation  (5) 

*The  concept  of  a complex  modulus  is  explained  in  reference  (5).  The 
addition  of  an  imaginary  term  to  the  modulus  allows  a representation 
that  includes  internal  damping  effects  for  pure  harmonic  motion. 
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/)  u 

c*"  = — , if  M is  real,  c is  the  velocity  of  sound  in  the  medium 
p = mass  density  of  material 
t * time 

o = complex  stress  in  the  medium 

For  a layer  whose  lateral  dimensions  are  large  compared  to  its  thick- 
ness (5): 


M - B + 

B = complex  bulk  modulus 
G = complex  rigidity  modulus 
For  (9),  we  assume  a solution: 


where 


u = UQ  c i(lJt  + 

Performing  the  appropriate  differentiations  on  (11),  and  substituting 
these  into  (9)  results  in  k = + £ , so  that  we  set: 


(11) 


u = U 


(12) 


Setting 

» =Sr  = Mr  ♦ • «» 

and  using  (10) 

„ = -Mf ju.«*f  -U+  c (14, 

Equations  (13)  and  (l4)  are  evaluated  at  the  input  and  output  sides 
of  the  plate:  f 


8 


*1  = ?! 


ttU 

V1  "at 


■ M gu 
ax 


NUSL  Tech  Memo 
No.  2133-648-69 


bc=0 


loie 


io)t 


..  10)  io>t 

-M  — e 
c 


[u-+  u+] 


[u--  ”♦] 


(15) 


(16) 


x=0 


*2  = ?2  = 


..  io>  io>t 
= -M  — e 
c 


[u  e1*-1  -u  e -^1 

L • c + c J 


X=1 


_ au 
v2  "at 


= iwe 


lojt  TT  nol  TT  „ -l  ojI 

U e U,  e — 

c + c 


(17) 


(18) 


k=i 


Equations  (17)  and  (l8)  can  be  solved  simultaneously  for  U e lwt  and 

io>t 


U+e 


, resulting  in: 

TT  . iwt  . 


1 — iual  f c I 

S’  e~ [V2-Tp2j 

IT  ^ _ 1 _ i^lL  C 1 

U+e  " 2U  e~[V2  +“m  P2J 

Substituting  (19)  and  (20)  into  (15)  and  (16): 


0)1 


iM  oil 


(UJ*  . . X1VX  (1^1.  1 

cos— ) p2  + (—  sin— ) 


ic  cal 


0)1  , 


V1  * M Sin  T > P2  + (cos  T] 


(19) 

(20) 

(21) 

(22) 
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Now: 

i?  ■#■£(£>■$  c*  ■ -C  = *P  <»i 

Substituting  (23)  into  (21)  and  (22),  and  writing  in  vector  form: 


pl 

I ■ 

10 1 

cos  

c 

iR  sin  — 

P c 

p2 

V1 

- 1 0)1 

— sin  — 
iR  c 

P 

0)1 

cos 

c 

v2 

Equation  (24)  is  a "typical"  four’  pole  representation  for  the  system 
of  Figure  7.  It  can  be  used  for  any  layer  in  a composite  system 
provided  that  M,  1,  and  p for  the  material  of  the  layer  are  sub- 
stituted into  the  four-pole  matrix.  A typical,  electrically  analogous 
T-network  is  obtained  below  by  using  equations  (5),  (6),  and  (7),  and 
is  shown  in  Figure  8. 
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Figure  8 

Cascaded  plates  are  represented  by  cascading  the  individual  T- 
networks  as  shown  in  Figure  9-  An  equivalent  single  T-network  can  be 
computed  by  using  Norton's  and  Thevinin's  theorems  directly  on  the 
electrical  analog  or  perhaps  more  easily  by  using  equation  (8)  to 
obtain  an  equivalent  single  four-pole  matrix,  and  computing  a T- 
network  from  this  using  equations  (5),  (b),  and  (7). 
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INRJT  AND  OUTPUT  NETWORKS 

The  electrical  analogy  of  Figure  9 is  not  complete  until  input 
and  output  networks  are  specified.  Figure  10  represents  a composite 
plate  upon  which  is  incident  a sound  pressure  pj. 


Figure  10 

There  is  a reflected  pressure  p_,,  and  a transmitted  pressure  p~ 
Comparing  Figures  10  and  7,  H 1 

Pi  = Pi  + Pr 
Pg  = Pi 

The  medium  on  the  input  side  of  the  plate  has  specific  acoustic  im- 
pedance R-^  = P-^  C]_  . The  specific  radiation  impedance  on  both  sides 
of  the  plate  is  assumed  equal  to  the  specific  acoustic  impedance, 
since  the  plate  is  of  large  lateral  extent.  The  specific  radiation 
impedance  on  the  output  side  of  the  plate  is  (6) 

^2  ~ ^ 2 ^2  ’ ^ VT 

is  the  velocity  associated  with  the  transmitted  wave  (6), 

p2  = P 2 C2  v2  = R2  v2 

Since  R2  is  real,  the  pressure-voltage  analog  of  the  output  is  as 
shown  in  Figure  11. 

Rt 


2 PIT  R1  V1 


Input  Network 


Output  Network 


Figure  11 
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Similarly,  denoting  by  vj  and 
incident  and  reflected  waves, 

PI  = 

% " 

Continuity  at  the  input  face  of  the  plate  requires  that: 

V1  = VI  + VR 

Hence: 

Pi  = Pi+Pr  = Pj-Ri  vR  = PI-R1(v1-vI)  = PI-R1  V-l+R-l  Vj-  = 2 Pj-^  Vj_ 

The  input  voltage  to  the  network  of  Figure  9 is  then  2 pj  - R^  Vj^  . 
This  can  be  realized  by  the  input  network  shown  in  Figure  11.  The 
total  network  for  the  plate  is  shown  in  Figure  12,  The  analogous 
quantities  of  interest  are  listed  in  Table  2. 


vr  the  velocities  associated  with  the 
respectively^ . 


R1  V1 


-R1  VR 


i — 'VW\ — o~ 
A 


2piO 


o 1 

Figure 

O 1 

12 

Acoustic 

Electrical 

Input  Velocity  Vj_ 

Current  through  R-^ 

Transmitted  Pressure  p^. 

Voltage  Drop  BC 

Reflected  Pressure  pt 

Voltage  Drop  AC  minus  \ 

source  voltage  current 

through  R2 

Output  Velocity  V2 

Current  through  R2 

Table  2 
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APPLICATION  OF  THE  FOUR-POLE  TECHNIQUE  TO  COMPUTE 
SOUND  TRANSMISSION 


Some  Specific  Examples 


Case  I - steel  plate,  thickness  small  compared  to  an  acoustic  wave- 
length (jjJl  < ^ ; 


M = ft  = ft  + 4/3  'o',  ( r,  = 0) 


Since  ™ is  small,  we  may  use  the  approximations 


tan  “i  = , esc  ^ = £_ 

2c  2c 


c wl 


Eef erring  to  Figure  8, 


Z-l  = Z2  = i Rp  = ipc  = iw£i 


2c 


2c 


pi  is  the  mass  per  unit  area  of  the  plate,  and  is  denoted  m . 

P 

Then 


Z-i  — Z0  = 


The  electrical  T-network 
13. 


2 -ift 

0)1 

^rorresponding  to  Figure  8 is  shown  in  Figure 


Z,  = -iP  -1  = “ i”  = _1  , 

J P wl  (j) 1 0)1  I 

l“(*> 


If  the  plate  thickness  is  on  the  order  of  1 inch, 

of  10"7,  while  nip  is  on  the  order  of  1.  Hence  we 

capacitor  of  Figure  13  and  the  analogous  circuit 
Figure  11  becomes  that  shown  in  Figure  14. 


2 


A » 


<0 

Figure  14 


C 

-a 
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1 , 

K is  on  the  order 
W 

may  disregard  the 
corresponding  to 


Two  important  acoustic  parameters  for  a plate  are  the  transmission 

oss  (TL),  and  reflection  loss  (RL).  These  are  computed  for  this 
case  2 


15 
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RL  = 10  log 


(i  + ±£.  ) 

\ <*>2  “p2  / 


Case  II  - Viscoelastic  Material, 

M = E + ^G  =^(1  + il) 

Using  the  same  approximations  as  Case  I, 


wl 

c 


< .1 


Zi  = Z2  = iw: 


2 


-ipc2  _ -iM  = + 1 


J3  wl 


w 


]_  wl 


In  this  case,  may  not  be  neglected,  and  the  analogous  T-network 
M 


Cases  I and  II  can  be  combined  to  yield  an  electrical  analog  for  a 
steel  plate  with  a thin  viscoelastic  damping  coating  on  the  incident 
side.  This  circuit  is  shown  in  Figure  16  where  the  damping  material 
constants  are  subscripted  with  d,  and  those  of  the  steel  plate  with  p. 
The  acoustic  impedance  of  both  sides  of  the  composite  is  assumed  to 
be  R. 
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Figure  16 


CONCLUSION 

Transmission  of  sound  through  layered  media  is  easily  analyzed 
by  the  use  of  the  four-pole  method.  The  basic  pole  matrix  (the 
matrix  of  equation  (24))  is  used  for  each  layer  of  the  composite 
structure,  by  substituting  the  appropriate  constants  into  aij. 
Formula  (8)  is  then  used  to  obtain  an  equivalent  four-pole  matrix 
for  the  entire  structure,  and  if  an  electrical  analog  is  desired, 
equations  (5),  (6),  and  (7)  are  used  to  obtain  the  T-network  that 
appears  in  the  analogous  circuit  of  Figure  11. 

Application  of  this  method  should  be  useful  for  analysis  of 
sonar  dome  construction,  baffle  design,  isolation  mountings,  and 
other  similar  problems. 


st  /// 

S.  M.  HU PI NS KI  ' 
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